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1. Introduction 

The strangeness and charmness contents of the nucleon are interesting aspects of nucleon 
structure. Despite the difficulty in measuring these quantities experimentally, the strangeness of 
nucleon still attracts much interest because it is believed to make a large contribution to the cross 
section between some dark matter candidates and ordinary matter [|]]. 

The strangeness and charmness of the nucleon can be defined as dimensionless fractions: 

•' mn ^ ■ ' 

n ^ m,{N\cc\N) .J 2) 



There are lattice QCD calculations with dynamical fermions by several groups [g[ ^ ^ |5|, |6||. Two 
different methods are used. One is to calculate (A'^l^^jA'^) = from the Feynman-Hellmann 
theorem. Another is to calculate the disconnected three-point functions 



, <N{t)qq{t')N{to) > -<N{t)N{to) X qq{t') > 

^ ' ' °^ < N{t)N{tQ) > 

where qq{t') and the nucleon sink N{t) are zero-momentum operators. We adopt the latter method 
to calculate the strangeness and charmness of the nucleon. To extract the form factor from the ratio 
of the three-point function and the two-point function, we use a summed ratio to enhance the signal 
and suppress excite-state contamination 

t-i 

R'it,to)= £ Rit',t,to) (1-4) 

f'=/o+l 



and do a linear fit for large t: 



R'{t,to)^^^^const. + t{N\ss\N) (1.5) 



2. Strategies for the computation 



2.1 Overlap valence fermion on domain- wall sea 



The overlap fermion action obeys chiral symmetry at finite lattice spacing and is free of 0{a) 
errors. It is shown that the effective quark propagator of the massive overlap fermion has the same 
form as that of the continuum [||]. The 0{m^a^) error, which is important in the charm region, is 
estimated to be small on quenched lattices |l^] and even smaller on the dynamical domain-wall 



sea [11] due to HYP smearing. Thus, it is shown that it can be used for both charm and light quarks 



on the 24^ x 64 lattice DWF configurations [12]. 

We use the valence overlap fermion on A'y = 2 + 1 domain wall dynamical configurations 
obtained from the RBC and UKQCD collaborations |[T^]. This is a mixed action approach. Since 



the valence is a chiral fermion, only one extra low-energy constant A^ix needs to be determined [14, 



15], which turns out to be small [12] 



2 



The Strangeness and Charmness of Nucleon From Overlap Fermions 



M. Gong 



2.2 The stochastic grid source 

We introduce Z3 and Z4 stochastic grid sources to gain more efficiency. A grid source is defined 



as 

77(f) = £0% (2.1) 

where ^ is a sparse regular grid of lattice sites on time slice f = 0, and 0, is the random phase on 
site /. The corresponding quark propagator is 

G{y,^)=D-\y,x)^{x) 

= £ diD-\y,x)5,j = £ eiG{y,T) (2.2) 

And the loop on site / is: 

L(/) = 0;G(7,T])^G(M) (2.3) 

2.3 The low mode substitution technique for the proton correlation functions 

The proton correlation function is 

Cpro,oniy,x;'^;G^"\G^''\G'^"'y) = e^^'e"'^'"' \tr{TG^'''^"''\y,x)G^'^^^^\y,x)G^"''^''\y,x 

+tr (^FG^")""' {y,x)^ tr (g^^^^^' {y,x)G'^"'^"=' (y,f)) J (2.4) 

where the G^"/'^^""' {y^x) is the uld quark propagator from xtoy with color indices a,a' and the G is 
defined [|l^] as {CGC^^Y with the charge conjugation operator C and C = C75 = 72/475. The trace 
and the transpose operations are defined on Dirac space. 

With the masses of the u and the d quarks set to the same value, G(^), G("), and G^"') are the 
same propagator. The correlation function C(Gi,G2,G3) is a linear function of each of its three 
arguments. 

With the Z3 grid source, the quark propagator is the sum of the propagators from different sites 
on the grid with different Z3 phases 

G = £ QiGi (2.5) 

and the proton correlation function is 

C{G,G,G) = C{Y,idiGi,Y,jdjGj,Y,i,dkGi,) 

= T,i,j,k ^ j ^kC {Gi,Gj, Gk ) 

llij,kSij,kC{Gi,Gj,Gk) 
^ LiC{Gi,Gi,Gi) (2.6) 

The limit is at that with infinite noise sources or infinite configurations. Otherwise, there is unbiased 
contamination by the propagators from different sites. 
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With deflation in the quark matrix inversion, we have the propagator with low-mode and high- 
mode parts as well: 

G= G" + G^ 

= G^ + L-ftGf (2.7) 



and the proton correlation function is 
C(G, G, G) = CiG" + L- diGf,G" + ejG'j,G" + ^ QkG, 



k, 

C{G",G",G") + ZiC{diGl diGf, diGf) 
+ ^iC{eiGf,G",G") +^iC{G",diGf,G^) +L-C(G^,G^, ) 
+ L-C(e,-Gf , diGf, G") + L-C(0;Gf , G", diGf) + L-C(G^, diGf, diGf) 
+ Li^jC{diGf, djG),G^) +Zi^jC{diGf,G^, djG)) + Zi^jC{G^, diGf, djG)) 

+ Li^7or j^kork^iC{diGf , djGp dkGf) (2.8) 

By design, the three quark propagators of the proton should start from the same site on the 
source grid. In this case, we can drop the terms with Y^i^j and Y.i^iorj^kork=^i which are all con- 
tamination without useful signal. For the pure high mode contribution and mixed G^ and G^, this 
is implemented by the Z3 noise. For the case with all low modes, we can implement it by hand, 
since we have the wavef unctions in G^. 

C{G,G,G) = CiG'i ,G'i ,G^) +L.C(0;Gf , 0,-Gf , 0,Gf ) 

+E/ G{diGf, G",G") + Zi G{G", diGf, G^) + D GiG^^ ,G" ,diGf) 
+ 1,- G(0,-Gf , diGf, G") + 1,- G(0,-Gf , G", diGf) + 1,- G{G" , diGf, diGf) 
GiG" + ZidiGf,G",G") 
+G(G^,I,- diGf, G") + G{G",G",Zi diGf) 
+ Li [GidiGf, diGf, G") + GidiGf, G",diGf) + G^G^" , diGf, diGf) 

+ G{diGf, diGf, diGf)] (2.9) 

With this low mode substitution technique (LMS) [|l^, [T8|], the proton correlation function 
acquires only errors from noises on the high modes and mixed G^ and G^. Similarly, one can use 
the low mode substitution technique on propagators from smeared sources. The excited states are 
significantly suppressed by the smearing technique and the signal-noise ratio is much better. 

2.4 The low mode average technique for the loop 

We use a similar technique to calculate the loop. To get more statistics, we use the low mode 
average technique, wherein we use the low mode loop on all the space-time sites instead of on just 
the grid sites. 

The loop for the scalar density is well saturated by low modes with light quark masses, but it 
is worse in the strange region and does not saturate at all in the charm region. 
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Figure 1: The comparison of proton effective masses. The fitting results for proton mass are 1 .13(14)GeV 
for point source, 1.08(5)GeV for LMS with grid noise source and 1.14(2)GeV for LMS on smeared grid 
source with folding. 




Number of eigenvectors 



Figure 2: The low mode contributions to the loop. 

3. Simulation details 

We use the overlap valence fermion on the 2 + 1 flavor domain wall fermion configurations 
from RBC. The lattice size is 24^ x 64 with \/a = 1.73(3) GeV and the sea masses are mi = 0.005 
and lUh = 0.04. We use 37 configurations for the calculation. 

We use a double grid source for the nucleon two-point functions. The sources are on the time 
slices f = and f = 32 with Z3 noises. The grid cell is (6,6,6) which means there are (24/6)^ =64 
grid sites on one time slice. 

We use partial dilution sources for the loop propagator. The time slices are diluted to 4 pieces 
and an extra even-odd dilution is adopted. Each diluted piece is a even-odd grid source with Z4 
noises. The grid cell is (4,4,4,2) on 4D space-time which means (24/4)^ x (64/2)/2 = 3456 grid 
sites are inverted simultaneously as sources. 

The strangeness with m„rf = 0.016 (which gives niji = 330MeV) for the nucleon and = 
0.067 for the strange quark shows /y^ = 0.048(15). The contribution from the low modes and the 
high modes of the loops are fj^^ = 0.041(12) and /r^^ = 0.003(5), which shows the strangeness is 
dominated by the low modes. 

The charmness with nic = 0.63 and mud = 0.016 shows fj^ = 0.029(43). The contribution 
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Figure 3: Left pane: the summed ratio for the strangeness content of nucleon. Right pane: the contribution 
from the low modes and the high modes of the loop to the summed ratio for the strangeness content of 
nucleon. 



1 p 

0.8 
0.6 

0.4 
0.2 

'> 
-0.2 
-0.4 - 




R'(i) 



* i i 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 



"III 



Low modefe 
High modei 



Figure 4: Left pane: the summed ratio for the charmness content of nucleon. Right pane: the contribution 
from the low modes and the high modes of the loop to the summed ratio for the charmness content of 
nucleon. 

from the low modes and the high modes of the loops are fjj^ = 0.042(12) and fjf^ = —0.016(45). 

Since we are using the multi-mass algorithm to calculate the propagators, the form factors for 
light quarks are calculated simultaneously. With both sea and valence rriud = 0.016, the discon- 
nected insertion contribution to the kN a term is 46(19) MeV. 

4. Conclusions 

The strangeness and the charmness contents of nucleon and the TiN o term are presented. 
The low mode substitution and the low mode average techniques are very helpful since the low 
modes play an important role for the scalar form factors. For the nucleon valence quark which 
corresponds to nij^ = 330 MeV, we find a three-sigma signal for the strangeness content of the 
nucleon, fj^ = 0.048(15), but less than a one-sigma signal for the charmness content where the 
culprit is the noisy high-frequency modes. 

This work is still preliminary with few configurations in one ensemble. We will proceed with 
more configurations on both 24^ x 64 and 32^ x 64 lattices with different sea masses. The variation 
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Figure 5: The summed ratio for the disconnected insertion contribution to nN G term, 
method can be adopted to reduce the error bars further. 
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